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ABSTRACT 

It is well known that the outer parts of QSO accretion disks are prone 
to selfgravity if heated solely by orbital dissipation. Such disks might 
be expected to form stars rather than accrete onto the black hole. 
The arguments leading to this conclusion are reviewed. Conversion of 
a part of the gas into high-mass stars or stellar-mass black holes, and 
the release of energy in these objects by fusion or accretion, may help 
to stabilize the remaining gas. If the disk extends beyond a parsec, 
however, more energy is probably required for stability than is avail- 
able by turning half the gas into high-mass stars. Small black holes 
are perhaps marginally viable energy sources, with important impli- 
cations (not pursued here) for the QSO spectral energy distribution, 
the metallicity of the gas, microlensing of QSO disks, and perhaps 
gravitational-wave searches. Other possible palliatives for selfgravity 
include accretion driven by nonviscous torques that allow near-sonic 
accretion speeds and hence lower surface densities for a given mass 
accretion rate. All such modes of accretion face major theoretical dif- 
ficulties, and in any case merely postpone selfgravity. Alternatively, 
thin disks may not exist beyond a thousand Schwarzshild radii or so 
(0.01 parsec), in which case QSOs must be fueled by gas with small 
specific angular momentum. 
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1 INTRODUCTION 

The outer parts of steady, geometrically thin, optically thick, viscously-driven accretion disks around QSOs 
are predicted to be selfgravitating: that is, the Toomre stability parameter Q [eq. (H)1 is much less than unity; 



equivalently, the midplane density is much greater than the Roche value ( Shlosman fc Begelman 1987 ). As 
shown in §^ for a standard a disk, Q falls to unity at a distance r ~ 1CP 2 pc from a central black hole of 
mass 1O 8 M0 accreting at the Eddington rate. The mass of the disk within this radius is much less than 
that of the black hole and hence able to fuel QSO activity for only a fraction of the black-hole growth time, 
^Edd ~ 5 x 10 7 Eq \ yr. This is one reason to suspect that QSO disks are larger than 10 -2 pc, but there is also 
indirect observational evidence. The spectral energy distribution of typical QSOs shows a strong peak in the 
rest-frame infrared, which is often interpreted as reprocessing of light from the central source by an extended 



warped disk (Sanders et al. 1989). Nearby Seyferts and other active galaxies, thought to be low-luminosity 
analogs of QSOs, sometimes have resolvable nuclear disks. Most spectacularly, VLBI observations of maser 
emission in NGC 4258 and NGC 1068 indicate disks on parsec scales and have been used to measure the 



mass of their central black holes (Nakai et al. 1993; Greenhill et al. 1995; Glreenhill & Gwinn 1997) 
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At least on larger galactic scales, it is believed that strong self-gravity leads to star formation (Martin & 



Kennicutt 2001). If that is true of QSO disks, there is a danger that most of the gas would form stars, leaving 



little to fuel the QSO. Several theoretical attempts have been made to modify the standard a-disk model so 
as to extend gaseous QSO disks self-consistently into the selfgravitating regime. The options include driving 
accretion with global bars or density waves so as to increase the radial velocity and lower the surface density 
needed to sustain a given mass accretion rate ( shlosman & Begelman 1989 ) ; heating or mechanically stirring 
the disk with embedded stars or black holes so as to raise the temperature and lower the density of the gas 
(Collin & Zahn 1999a 3); or allowing the disk to fragment into colliding gas clumps whose epicyclic motions 
stabilize it against selfgravity, at least on scales larger than the clumps themselves (Kumar 199£). 

Our purpose is to re-examine selfgravity in QSO accretion disks with an emphasis on dynamical and 
energetic constraints. These constraints are most severe for massive and luminous systems, so our interest 
is in black-hole masses M > 1O 8 M0 and luminosities close to the Eddington limit. 

The plan of the paper is as follows. §2 reviews selfgravity in steady a disks, including irradiation from 
the central source (§ [2.2]), and enhancements to a by local gravitationally-driven turbulence in a moderately 
selfgravitating disk (§2.3). §3 considers disks that are stabilized by additional heating beyond that due to the 
dissipation of orbital energy, but in which angular momentum continues to be transported by an a viscosity. 
The effective temperature of such disks falls more slowly than the standard relation T c s oc r -3 / 4 ; in contrast 
to the usual situation, the total energy released per unit mass accreted is strongly dependent on the outer 
radius of the disk (^). §4 briefly explores various alternatives to viscous thin-disk accretion. Schemes that 
involve thin disk but invoke faster-than-viscous angular-momentum transport include global spiral waves 
(§4.5) and magnetized disk winds (§4.3). In every such case, we argue that the accretion velocity is bounded 
by the sound speed. Other alternatives explored in §4 are less like thin disks: quasi-spherical flows (§ [4.l[ ), 
collisional star clusters (§4.2) and clumpy disks ( §|4.6| ). 

All of these alternatives to the standard alpha disk face severe theoretical difficulties, or else seem 
unlikely to permit a centrifugally supported accretion flow beyond ~ 0.1 pc. We conclude in §5 that the 
gas probably does not circularize beyond this radius and must be supplied to the nucleus with low specific 
angular momentum. 



2 STEADY ALPHA DISKS 



For the time being, we assume steady accretion and neglect winds, so that the mass accretion rate is constant 
with radius, and advection of angular momentum is balanced by viscous torque. 
1. The mass accretion rate and viscosity parameter are related by 

M = inaffcln' 1 ^, (1) 

in which c s = \fpfp is the isothermal sound speed at the disk midplane, f2 = (GM/r 3 ) 1 ^ 2 is the orbital 
angular velocity, E is the surface mass density, and j3 = p gas /p is the ratio of gas pressure to total pressure at 



the midplane. The dimensionless Bhakura & Sunyaev (1973) viscosity parameter a is not necessarily constant 
but presumably < 1. The mechanism likely responsible for the viscosity of most disks is magnetorotationally- 
driven turbulence, for which simulations indicate a — 10~ 3 — 10 _1 (Balbus & Hawley 1998). As discussed 
below, where the disk is selfgravitating, a may be as large as ~ 0.3. Therefore, we often write ao.oi = 10 2 a 
or ao.3 = a/0.3. The parameter & is a switch that determines whether the viscosity is proportional to gas 
pressure (6 = 1) or total pressure (6 = 0). It is well known that radiation-p ressure-dominated regions of a 
disks have viscous instabilities in the latter case (Lightman & Eardley 1974), but it is possible that 6 = in 



an average sense. Although the average surface density E would be lower and the corresponding Q higher 
for 6 = than for 6=1, the viscous instability is expected to produce overdense rings in which these trends 
may well be reversed. 

2. If the disk were heated by viscous dissipation only, then its surface temperature would be 

T eff = r^^v 74 (2) 



V87tcr 
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2.9 x lO 3 (M 8 M ) 1/4 (r/l(T 2 pc)- 3/4 K 



w 6.2 x 10 5 ( I I — I K. 

\£ .iM 8 J \R S J 

We have introduced the abbreviations Mg = M/(1O 8 M0), Mq = M/(1Mq yr _1 ), dimensionless luminosity 
l E = L/Le, radiative efficiency s = L/Mc 2 = O.leo.i, and Schwarzshild radius Rs = 2GM/c 2 ~ 10 _5 Ms pc. 
The mass accretion rate can then be written as 

M = h « 2.2 eo.l h M S M , (3) 

where k CiS . « 0.4 cm 2 g -1 is the electron-scattering opacity. 

3. If viscous dissipation occurs mostly near the midplane, and vertical transport of heat is by radiative 
diffusion, then the midplane and surface temperatures are related approximately by 

r 4 -fr V (4) 

Simulations of magnetorotational turbulen ce indicate that the v ertical scale height of the turbulent dissipa- 
tion is larger than that of the gas density ( Miller fc Stone 200C ) . Hence eq. (^) may somewhat overestimate 
the midplane temperature, in which case selfgravity may extend to even smaller radii than estimated below. 



1/4 / \ -3/4 



l E V 7 * r 



2.1 Onset of selfgravity 

The Toomre stability parameter for keplerian rotation is 



(5) 



^ TtGE 27tGp ' 

and local gravitational instability occurs where Q < 1. We have taken E = 2hp with h — c s /fi, which is 
approximately the disk half-thickness when Q > 1. Combining eq. (JHJ) with eq. (|l|), one obtains the simple 
relation 

GMQ = 3a/3 b c 3 B . (6) 

For future reference, note that in a flat rotation curve, u c i rc = rQ. — constant, the numerical factor 3 is 
replaced by 2\/2 : a distinction too fine to matter. 

To see why selfgravity is inevitable, consider the radiation-pressure-dominated case, f3 <C 1, so that the 
isothermal sound speed c 2 = 4aT 4 /3cp. From eqs. (|) & (§), T 4 = 3reEfi 2 M /lQna. Eliminating T between 
these relations and using E = 2c s p/Q leads to 

c s = ^ = ^— flRs (P < 1), (7) 

whence h — (Z E /e)-Rs =constant for k = k c . s .. Using this for c s and eq. (Bh for M in eq. (o) leads to 



3a/3" /U 2 , 3 ( k c . s c a \ [Rs\ m 
8nV2 [ej K { G 2 M )\r ) 



(/3 « 1), (8) 



in which k = k/k . s .. Hence disks around more massive black holes are more prone to selfgravity. For b — 0, 
there is a characteristic mass above which an Eddington-limited disk would be selfgravitating even at its 
inner edge. This mass is enormous (~ 10 19 Mq for ao.oi = eo.i = 1), but even for realistic black-hole masses, 
Q < 1 at radii greater than 

Notice that a has been scaled to 0.3 rather than 10~ 2 because of the expected enhancement by gravitational 
turbulence (see below). 
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We now examine the assumption that <C 1. It follows from eqs. (^), and (^) that 

al-(b/2) /o \ 3/4 / \ 1/2 /, 4 9 2 \ I/ 4 / \ -3/4 

g 1 7 J _ 5/2 /3\ / ea\ n -3/4 I k B c V. s . \ ( r \ ' 



d-^ - w Q T^wJ U> M (10) 



0.39 [£H^£ q-3/4 (t7t ^) m - 



V2 , / r \ -3/4 

^lOOife. 

Eq. ([To]) does not depend upon eqs. (^) & (Q): that is, it does not assume that radiative losses are balance 
by viscous dissipation. (The opacity k . s . enters the formula only because it is used in the definition of the 
Eddington ratio l E ). Hence eq. ([To]) remains true even when auxiliary sources of heating are present in the 
disk. For viscously heated, radiatively cooled disks, /3 is given by eq. (A3) in the Appendix. 



2.2 Irradiation 

Although flared or warped outer regions can be warmed somewhat by irradiation from the inner parts (where 
most of the total disk luminosity originates), we now show this effect is not enough to stabilize the disk at 
r ~ r s . g .. 

To do so, irradiation must raise the midplane temperature T, not just the surface temperature T e ff, 
substantially. As a result, the disk will be nearly isothermal from surface to midplane. The vertical pressure 
gradient, and therefore the disk thickness, will be dominated by the gas even if /3 <C 1: that is, h 2 « p gas /pfi 2 . 
If the effective viscosity derives from local magnetorotational turbulence, then it probably scales more directly 
with thickness than with sound speed since MRI instabilities do not require compression. Then in all relevant 
respects the disk behaves as if /3 = 1 and c 2 s — > p gas /p = k^T/m. From eq. (^|), it then follows that the 
minimum temperature that the irradiation must provide to ensure gravitational stability is 

T =, =1 = f f^^^.exlO^^) 273 K. (11) 
\ 3a / \ Q o.3£o.i / 

(To be conservative, we have scaled to a = 0.3; for a = 0.01, this temperature would be an order of 
magnitude larger). On the other hand, the temperature of the disk in equilibrium with radiation from the 
central source is 

where 6 is the angle between the local normal to the disk and the radial rays. For a flat disk, cos (9 ~ 
max(/i, Rs)/r at r > _Rs, but even for a severely warped disk with cos# ~ 1, T cq <g; Tg =l g = i at r/7is 3> 
44(ao.3eo.i) 4/,3 ^E 5 ^ 6 Afg~ 11 ^ 6 - Since this last inequality certainly holds at all r ^ ?"s.g., we conclude that 
irradiation is not important for the selfgravity of the disk. 



2.3 Local gravitational turbulence 



It has occasionally been proposed that a partially selfgravitating disk may transport angular momentum by 
spiral waves (e.g. Cameron 1978| ; Paczyhski 197ij ). Others have suggested that gravitational instabilities are 
intrinsically global and therefore not reducible to a local viscosity (Balbus & Papaloizou 1999). It may be 



that both opinions are sometimes correct, depending upon the ratio of disk thickness to radius, since the 
most unstable wavelength in a Q = 1 disk is ~ h. We consider local gravitational turbulence here, since it 
can be accommodated within the a model; a (much larger) upper bound on angular momentum transport 



by global spirals is given in §4.5 



By careful two-dimensional simulations, Gammie (2001) finds that gravitationally-driven turbulence can 
be local and can support a approaching unity; in fact, in the absence of any other viscosity mechanism, his 
disks self-regulate themselves so that 

(13) 



(72D - l)fitth ; 



© 2002 RAS, MNRAS 000, 



Selfgravity and QSO disks 5 



where t th = ^k B T/aT* B is the local thermal timescale, and 72D is the 2D adiabatic index: 



72D 







dlogE 



Ios 



) dz 



[Our eq. (|l^) differs from Gammie by a factor 72D because we define a in terms of isothermal rather than 
adiabatic sound speed.] Gammie finds that self-regulation fails and the disk fragments if Q.tth < 0.3 for 
72D = 2. The latter result is not entirely surprising: in the absence of pressure support, the natural timescale 



for collapse is t^ y 



Q , and one would not expect collapse to be prevented by thermal energy unless 



tth > idyn- The converse statement, that fragmentation can be indefinitely postponed if Qt t h > 0.3, is not 
at all obvious but is supported by Gammie's simulations. 

We are not aware of any direct numerical simulations of disks that are unstable to both gravitational 
and magnetorotational modes at the same radii, as is likely to be the case for QSO disks (Menou & Quataert 
2001). It would be interesting to explore this, as there might be a synergy between the two instabilities. 

QSO disks are expected to be very thin in the regions of interest to us, so that Gammie's results may be 
applicable. Eq. (Q) implies h = l E Rs/6e in a radiation-pressure dominated disk, and hence h/r < 0.003 at 
the innermost radius where Q ~ 1 [eq. (g)]. We expect that from this radius outwards, a will rise smoothly 
from the value supported by magnetorotational turbulence (perhaps a m .h.d. ~ 10 -2 ) up to the maximum 
allowed by gravitational turbulence, a gra v,max ~ 0.3, in such a way that Q ^constant. It follows from eq. (^) 
that the ratio of outer to inner radii of this region is rather modest: ~ (a g rav,max/o!m.h.d.) 2 ' /S ~ 2. At still 
larger radii, additional sources of energy are required in order to prevent catastrophic fragmentation. 



3 CONSTANT-Q DISKS 



In this section, we postulate that some heating process arises that maintains the disk on threshold of 
gravitational instability, i.e. at a constant Q « 1, and work out some consequences for the structure of the 
disk. In 



3.1 



we discuss the radial dependence of midplane temperature, density, and (3 [already given by 



eq. (jlOp] . In §3.2, we estimate the minimal amount of energy in excess of viscous dissipation that must be 
supplied to maintain constant Q. 

For our immediate purposes, the details of the heating process are not important as long as it provides 
a stable feedback that regulates Q. Angular momentum is assumed still to be viscous, though perhaps with 
an enhanced a in eq. ([j]). The effective temperature of the disk may no longer be determined by equation 
(S), however, because viscous dissipation may not dominate the thermal energy budget. 

Although they may seem artificial, these assumptions are probably appropriate for the disks of spiral 
galaxies, and for the local interstellar medium in particular. The local Galactic magnetic field is consis- 
tent with simulations of magnetorotational turbulence: viz. a somewhat sub- equipartition strength, a pre- 
dominantly toroidal orientation, and fluctuations comparable to the mean ( Brandenburg et al. 1995 ). It 
is plausible therefore that there is a nonzero average magnetic stress (B r Bg) /47t = ap gas that systemati- 



0.3m_H cm 



cally transfers angu lar momentu m outwards ( Scllwood 

3 (gpitzcr 1978), and circular velocity Vb 



Balbus 1999[ ). Ta king p gas /fc B g 2000 cm~ j K 



200 km s _i ( [Binney fc Merrifield 1998| ), the 
implied radial drift velocity v r w — ap ES , B /pVo ~ —0.3 a km s -1 is small enough to have escaped detection. 
Perhaps coincidentally, eq. ([]) predicts M ~ 3 x 10~ 2 ao.iQ -1 Af© yr _1 , about half the Eddington rate for the 
Galaxy's 2.5 x 1O 6 M0 central black hole. The implied viscous heating rate ap gas £2o ~ 2 x 10~ 29 Q0.1 erg cm -3 
is negligible compared to the inferred radiative cooling rate of the gas, ~2x 10 _26 ergcm -3 (Spitzer 1978). 
Presumably, the temperature of the ISM is maintained by stars. An important difference between the the 
local ISM and QSO disks is that the former is very optically thin, especially to absorption, which means 
that the energy input from stars is inefficiently radiated. 



3.1 Density and temperature 

The midplane density in a constant-Q disk follows from eq. Q: 
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?~ = ^Q-i(^)\™->, (14) 

so that the density at r s . g . [eq. (§)] is ~ 10" 8 M 8 " 4/3 gem" 3 . The ratio (3/(1 - (3) = 4aT 3 /Sep is determined 
by eq. (|lo|). The temperature itself is (for /3 <C 1) 

1/12 / \ 1/6 / R „ x 3/4 

T = 2~"" 



-7/6 / 3\ / l E \ 19/12^-5/12 -1/4 -1/6 Jf-1/3 f ^.V 



1/6 / p \ 3/4 

M -l/3 £tS K iffc = . ( 15 ) 



' E ^| M" 1/3 ( — s 
Mo.iQ 1 / 2 J 8 I r 



t = 2 - i / 3 (|) i/6 ^^^y / % 5 / 6 G i /v- i /v,- i/3 (w^- i/3 1 J 

« 1.4 x 10 6 f M^y 7 VM 1/2 K if 6=1. (16) 
The surface density is 



E = 



/ E c 8 \ 1/3 3/2 



12n 2 V2 eaQ 2 G 4 K e . s M 2 J \ r 



I \ 1/3 , /Ro\ 3/2 



, , M 8 " z/J ( — I g cm" 2 if b = 0; (17) 
<>j.3£o.iV 2 / V r 

2 / 3 / Q 7 \ 1/6 



o / o; \ 2 /3 / o 7 \ I/ / d 

2 / 2/ E m \ I 6<JC \ /its 



3 V £ a feces. / \7xQG) v r 



5.4 x 10 8 Q- 1/6 (^^) 2/ V^') gem- 2 if 6 = 1. ( i,s ) 



ao.3£o.i 

For the nominal values of the parameters shown, E would fall to lg cm -2 at 40 pc for 6 = and 5kpc for 
6=1; however, in the latter case h/r rises more rapidly with r and is already ~ 0.5 at r — lpc. 

Beyond 10 4 — 10 5 i?s ~ 0.1 — lMspc, the above formulae predict T < 5000 K, so that the opacity 
K <C Kc.s. and the disk becomes optically thin. (This assumes Mg = l E = eo.i = c*o.3 = 1- Dust will raise 
opacity again at T < 1700 K.) The disk must then be supported by gas pressure, notwithstanding eq. ( |i"c| ) 
which presumes that the radiation is trapped. But at ft = 1, the minimum temperature for gravitational 
stability is ([n]), which is about an order of magnitude larger than predicted by the formulae above. So in a 
marginally gravitationally stable disk, there must be an extended region where the temperature adjusts itself 
within a limited range (5000 K < T < 10 4 K) so that the disk is marginally optically thin. At the low densities 



relevant here, the maximum opacity is K ma x ~ 10k CiS . and is achieved at T ~ 10 4 K (Kurucz 1992; Keady & 
Kilcrease 200C| ). Hence the outer edge of the region in question should end at E « « 0.3 cm g 1 which 



occurs (assuming Q = 1 and T = 10 K) at 
rthin ~ « 170M 8 pc, (19) 

provided of course that the disk extends to such large radii. We have taken the circular velocity as <7buige\/2 oc 
M [see eqs. (|5|)&@] rather than ^GM/r, since r thin lies outside the black hole's sphere of influence 
[eq. (§)]. 

The parts of the disk beyond r t hi n would radiate predominantly in optical emission lines with velocity 
widths ~ cTbuigc, so that they might be identified with the QSO narrow-line region. However, we will now 
see that the energy required to maintain the disk at constant M and Q > 1 all the way out to r t hi n would 
be prohibitive. 
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3.2 Energetics 

We define a local disk efficiency by 

For a viscously heated disk in newtonian gravity, e c 2 reduces to the local binding energy per unit mass, 
GM/2r = c 2 R s /4r [see eq. (Eh], which is largest at small radii; given a torque- free inner boundary at 
t — r m i n , the global efficiency is s = e'(r m i n ). But the constant-Q disks require additional energy inputs, so 
that e' is generally larger than Rs/4r and tends in fact to increase with radius. 



Using the results of §3.1 for the midplane temperature T, and assuming that vertical radiative transport 



obeys eq. we find that 

'm _ 711/3 ( Qe2 GRs V /3 (i V /a 

£[r) 2V631/3 K ^ a p C 2 Kc s ) y Rs ) 

=-2 \ 1/3 / „ X 1/2 



ft-iM^ 8 ( \ if 6 = 0; (21) 



kIL 3 Q^e 1/3 (G\ 5/6 fm\ 2/3 



1/6 



2V635/6 K a 2/3;V3 V c / V fc 

1.1 x 1(T 2 ( Q5/ 2 2£ , 01 ) «" 1m b f —^r ) if 6 = 1. (22) 



n 



2 , 

0.3'E 



; E J V io 5 i? s 



It would appear from these relations that e' — * 00 as n — * 0, but they are not valid when the disk is optically 
thin. If the absorption optical depth r = < 1, then T^j w rT 4 rather than T 4 /t. So in fact e' — > as 
k — > 0. Hence the largest radius at which eqs. ( |2~l| ) or ( pj[ ) is valid is the smallest among the radii r rt hi n , Rn, 
and r ou t) where the first is the radius at which a Q = 1 disk would become optically thin [eq. ([l9])], the 
second is the radius within which the black hole dominates the potential, 



i?N = — 2 « (6 ± 1) A4 U * ±U PC , (23) 

/0 bulge 



and the last is the actual outer edge of disk, which depends upon the angular momentum of the gas supplied 
to it . The final expression in eq . (|23| ) has used the recently reported M — <7bui gc relation ( jGcbhardt et al ' 



200C 



Merritt & Ferrarese 2001 



see §1.2 below). Actually, eq. (E3j) probably underestimates the distance to 



which the black hole dominates the circular velocity, because the density profile in the cusp of bright bulges 
is considerably shallower than r~ 2 . 

At r = 10 pc « 10 6 Mg 1 R s , equations @ and (§2|) imply s' w 0.003 Af 8 " 1/6 and e' w 0.1, respectively. 
The former efficiency is barely compatible with thermonuclear burning even if all of the disk gas is processed 
through high-mass stars. The latter is unsustainable by stars and in fact comparable to the efficiency of the 
central engine. 



4 ALTERNATIVES TO THIN DISKS 

As shown in §|^, a geometrically thin, optically thick accretion disk in a typical high-luminosity QSO would 
be selfgravitating at radii r > 10 2 - 10 3 i?s, or 10~ 3 - 10~ 2 pc, where R s = 2GM/c 2 is the Schwarzschild 
radius of the central black hole. Large global radiative efficiency, e > 0.1, probably requires a thin disk near 
Rs, but it does not much depend upon the nature of the flow at large radius. So in this section, we consider 
whether the selfgravitating part of the disk can be replaced by some other form of accretion. 
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4.1 Hot, quasi-spherical flows 

These have been proposed as models for accretion at very low accretion rate (M) and low radiative efficiency 
e = L/Mc 2 (Rees et al. 1982; Ichimaru 1977; Narayan & Yi 1994). At sufficiently low density, the gas 



cooling time and ion-electron thermal equilibration time are longer than the accretion time, so that the 
ion temperature (Ti) is approximately virial and the thickness of the disk is comparable to its radius. One 
supposes that angular-momentum transport is efficient (viscosity parameter a ~ 1), or that the angular 
momentum is very small to begin with, so that the accretion velocity (v r ) is comparable to the free-fall 
velocity. The large v r and Ti combine with the low M to make the density low, as required. It is unclear 
just how low M must be for this mode of accretion to sustain itself, because angular-momentum transport 
is not well understood and collisionless processes may enhance the thermal coupling of ions and electrons. 

A quasi-spherical flow is not viable when the luminosity is close to the Eddington limit, however, 
because of inverse-Compton cooling. For spherical free fall onto a source of luminosity L — /s^Edd, the 
inverse-Compton cooling rate of free electrons is 

1/2 



2otL 
3nr 2 m e c 



_ 8 rn^ I Rs 
9 m e V r 



We have introduced the free-fall time of a radial parabolic orbit from radius r, 



ts 



3 V 2GM 



1/2 



-V2 r 3/2 
1 8 pc a i 



(24) 



(25) 



in which r pc is the distance from the source in parsecs and Ms = M/(10 8 Mq). Actually, radiation pressure 
increases the free-fall time by (1 — Z E ) -1 / 2 , but this factor is neglected for simplicity. From eq. (^), tc < iff 
at r < 25/|Ms pc. Hence the electrons assume the color temperature of the central source (or even less, see 
below): Tc < 10 s Mg -1 '' 4 . This is much less than the virial temperature, T vl 
interest to the present paper. 1 The electron density is 



10 7 M$r P c , at all radii of 



£ I 3C(TT iff ' 



(26) 



so that the electron- ion equilibration time due to Coulomb collisions alone (Spitzer 1978) is much shorter 
than the flow time: 

3/2 



^ « 2 x Kr 5 d, 
iff 



f T e 



1 10 5 K 



Of course T e ~ 10 5 K is the peak of the cooling curve (Spitzer 197S ) . If n e obeys eq. (|2q), then radiative 
cooling is actually faster than inverse-Compton cooling: 



tc V r 



trad 



(T e = 10° K). 



In short, both the ions and the electrons of a quasi-spherical flow onto a near-Eddington QSO cool in 
much less than a free-fall time. A thin disk will form unless the specific angular momentum of the flow is 
negligible. 



4.2 Collisional stellar cluster 

Dense stellar clusters have occasionally been nominated as precursors to QSO black holes, either by relativis- 



tic collapse (Zel'Dovich & Podurets 196f; Shapiro & Teukolsky 1985; Ebisuzaki et al. 2001) or by collisions 



among non-degenerate stars (Spitzer & Saslaw 1966; Rees 1978). Our interest, however, is in stellar collisions 



as the main source of fuel for an already very massive black hole. This has been studied by 



McMillan et al 



P~[The timescale for Poynting-Robertson drag is (m p /m e )tc, hence ^> tg at r 2> i?s> *- e -- 
not remove angular momentum from the gas fast enough to prevent it from circularizing. 



the radiation field does 
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(1981) and [Uarionov & Romanova (1988), among others, who show that in order to supply a 10 s Mq black 
hole at its Eddington rate, the velocity dispersion of such a cluster must be > 10 3 km s _1 . To demonstrate 
the robustness of this conclusion, we will make some oversimplified but conservative estimates here. 

In order to provide high radiative efficiency (e), stars must be disrupted and their gaseous debris 
circularized before accretion. If Mb.h. 10 s A/q, main-sequence stars scattered onto loss-cone orbits are 



likely to be tidally disrupted rather than swallowed whole (Hills 1975). About half of the tidal debris is 



unbound and promptly escapes from the black hole (Lacy et al. 1982; Evans & Kochanek 1989), and much 



of the remainder is likely to be swallowed at low radiative efficiency before the gas circularizes (Cannizzo 



st al. 1990 



Ayal et al. 200C). We limit our discussion to Mb.h. > 1O 8 M0, as required for the most luminous 



QSOs, and assume that stars are disrupted by stellar collisions. We ignore loss-cone effects because stars 
swallowed whole do not contribute to the QSO luminosity. 

The total collision rate among N stars forming a cluster of structural length a is 



N 



Ca 7 N 2 {GM)~ 3 Rl, 



(27) 



where C is a dimensionless coefficient, a the root-mean-square velocity dispersion in one dimension averaged 
over the cluster, N the number of stars, and R, ~ Rq the radius of an individual star. M is the total mass 
that determines a via the virial theorem, so that if m* is the mass of individual stars and M* = Nm* then 
M ~ Mbh + ^M, , the factor of 1/2 being needed to avoid double-counting the gravitational interactions 
among the stars. We take 4nR 2 for the collision cross section; this allows for grazing collisions that probably 
would not disrupt the stars (Spitzer & Saslaw 1966), but it will lead to a conservative estimate of a. 



Gravitational focusing is unimportant at the high velocity dispersions relevant here. 

A second relation for N follows by requiring the collisional debris to sustain the QSO at a fraction 
l E = L/L-e of its Eddington luminosity: 



2M,N 



1-e 4nGMm p 

£ C(TT 



(28) 



Eliminating N between eqs. (j27j) and p8[ ) leads to 

_ 7 _ 2nl E G 4 m p m„ / M b . h .M 3 " 
° ~ Ce a T c ~Rj \ Mi 

The term in parentheses achieves its minimum value 27M b ! . h ./16 at M* = 4M b . h .- 

The coefficient C depends upon the density profile p(r) of the stellar cluster. It can be arbitrarily large 
if p(r) is sufficiently steep, but then the collision rate is dominated by stars at small radius. These tightly- 
bound stars represent only a small fraction of the cluster mass and would be consumed in much less than 
the growth time of the black hole, unless their total mass (Mt.b.) is > Mb.h., in which case the tightly- 
bound population is substantially self-gravitating and we redefine M* = M t .b.- One might suppose that 
M t .b. Mb.h. and that the tightly-bound stars were continuously replenished by two-body relaxation from 
a reservoir of more weakly bound stars; but the larger cluster would then have to expand to conserve energy, 
reducing the collision rate and the fueling of the QSO. For these reasons, we consider clusters for which 
C is dominated by stars near the half-mass radius. For example, if the stars have isotropically distributed 
orbits with a common semimajor axis a in a keplerian potential, then p(r) oc \J (2a/r) — 1, a 2 — GM/3a, 
and C w 196. This leads to 



1/7 



1/7 



Ml 11 km s" 1 , 



(29) 



where gr» is the stellar surface gravity. Because of the one-seventh root, the result is not very sensitive to our 
assumptions. A Plummer sphere of the same total mass, M* = 5 x 10 8 Mq, yields a ~ 730 km s -1 . 

Equation (^9|) can be compared with recently-discovered empirical relations between inactive black 
holes — presumably QSO relics — and their host bulges. Gebhardt et al. (2000| ) find 



Mb 



1.2(±0.2) x 10 8 



200 km s" 



Mr: 



(30) 
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where <r is the line-of-sight velocity dispersion at one effective radius. Merritt & Ferrarese (200f ) use the 
central velocity dispersion, which is usually little different from a c : 

M b . h . = f .30(±0.36) x 10* M . (31) 

The scaling exponent dlogM/dloga = 3.5 implied by eq. ( p9| ) is similar to the empirical ones (j?o|) & 
([jl|), but the normalization is very different. Extrapolated to 700 km s _1 , the empirical relations predict 
Mb.h. > 10 10 Af Q instead of 1O 8 M . 

Therefore, if QSOs were fueled by dense stellar clusters, these clusters must have been an order of 
magnitude more tightly bound than the surrounding bulge, and they would have been a dynamically distinct 
stellar component. There seems to be very little trace of this tightly-bound stellar population in present- 
day bulges. How would such a component form? A likely possibility is gaseous dissipation followed by star 
formation. As will be seen, accretion in a thin viscous disk may lead to just such a result. We have discussed 
fueling the QSO by a disk and fueling it by stellar collisions as though these were mutually exclusive 
possibilities, but perhaps the two occur in concert. 



4.3 Wind-driven disks 

In principle at least, a magnetized wind can remove angular momentum from a thin disk rather efficiently. 
Compared to a viscous disk of the same sound-speed (c s ) and accretion rate (M), a wind-driven disk might 
have an accretion velocity that is larger by a factor ~ rQ./c s = r/h. The surface density would be corre- 
spondingly reduced, as would the tendency toward self-gravity. 

If viscous transport can be neglected, then under steady conditions, 



M aC cfir 2 = -y BzB^r^dr'. (32) 

r min 

The lefthand side is the angular-momentum flux carried inward by the accreting gas. The mass flux through 
the disk M acc , is now a function of r even in a steady disk because of the mass loss through the wind, 
Mwind -C Macc The righthand side of eq. ( [32]) is the rate of loss of angular momentum by maxwell stresses 
exerted on both faces of the disk, assuming that B$ is odd in z while B z is even. The overbars denote an 
average over fluctuations in azimuth and time. The minus sign appears because we define M acc > for 
inflow. The corresponding equation for a viscous disk is 



h 



M acc f7r 2 = -ir 2 j B r B$, dz. (33) 

-h 

The magnetic part of the angular-momentum flux emerges through surface areas 4n(r 2 — r 2 ^) and Anrh in 
eqs. J32] ) & (|33|), respectively. Since the former area is larger than the latter by a factor r/h, the field needed 
to drive a given M acc is smaller by h/r in a wind-driven disk than in a viscous disk. 

If the effective viscosity of an accretion disk is magnetic, then the energy density of the field is at most 
in equipartition with the gas, viz /3 ma g = 87tp gas /B 2 > 1, because a super-equipartion field would shut off 



the crucial magnetorotational instability (Balbus & Hawley 1995). We now give two arguments to show that 



even in a wind-driven accretion disk, the field must also be at or below equipartion. 
Vertical hydrostatic equilibrium requires 



E?2 1 d2 n2 

8n +Psas 



pQ 2 z dz . 



We assume that the field has dipolar symmetry, so that B z is approximately constant with z on the scale 
h <C r, whereas the horizontal components vanish at the midplane. Since the righthand side above is positive, 
and since p gas (z = h) <C Pgas(z = 0), it follows that 
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Bl + B% 

871 



(* = ±ft) <Pgas(z = 0). 



(34) 



Henceforth _B r and are evaluated at 2 = ft, and p gas at 2 = 0. As pointed out by Blandford & Payne 



(1982), in order that centrifugal force should drive the wind outward, the poloidal field lines must make an 



angle of at most 60° with the surface of the disk, so that B r /B z ^ %/3- Therefore 

B 2 ^ 4B 2 +B1 

8^ 55 87T ^ 4Psas ' 

The magnetic energy density is less important to us, however, than the azimuthal force per unit area on the 
disk: 



B z B<f, 



47T 



B r B^, 



in 



^ V ^Bl + B ±< ^ 



8tt 



< V3p gas . 



Using the inequality in eq. (p2|), together with M a , 



-2nY,v r and p gas ~ T,c 2 /2h, we have 



fir 3 £ 



y> 2 r 2 

z -' c s 1 1 \ 12,1 AT L s 
_( r ) r dr 



(35) 



(36) 



So the inflow could be marginally supersonic. For a viscous disk, on the other hand, \v r \ w ac 2 /f2r, which is 
strongly subsonic as long asa<r/ft. 

The angular momentum extracted from the face of the disk by the field must be carried off by the wind. 



It is problematic whether a strong wind can be launched from the disk (Ogilvie & Livio 2001), and whether 
rapid wind-driven accretion is stable (Cao & Spruit 2001). Apart from these difficulties, consideration of the 
magnetic flux, 



$(r) 



B z (r',0) 2nr'dr'. 



leads to important constraints on the field strength and accretion rate. Presumably $(r) should not change 
secularly. If B z is predominantly of one sign, then the increase of |<E>| by advection must be balanced by 
diffusion of the lines through the inflowing gas. The drift velocity of the lines Vdrift = ~ v r ~ T] e «/r if 
B z varies on scales ~ r, where r) e s is the effective diffusivity of the gas. In significantly ionized disks, 
the microscopic diffusivity is negligible, so ?7 e ff is due to turbulence, and one expects r] e « = a mag c s /i with 
a ma g 1 for much the same reasons that the effective viscosity v c g — a ma gC s /i with a v i sc < 1. Hence for 
a steady wind-driven disk threaded by net magnetic flux, \v r \ < (a mag /i/r)c s , which is probably very much 
less than the upper limit (j3(j) and comparable to the accretion velocity of a viscous disk. Alternatively, 
the net flux could be essentially zero if B z changes sign on scales r. In the latter case, the higher flow 
speed ( |36[ ) may be achievable. But so irregular a field would probably have to be sustained by dynamo action 
within the disk rather than inherited from whatever region supplies the accreting gas. This probably requires 
magnetorotational instability (henceforth MRI) and gives another argument for a sub-equipartion field. 

To summarize this subsection, accretion driven by magnetized winds is even less well understood than 
viscous accretion but might allow substantially higher accretion velocities and lower surface densities, perhaps 
by factors up to ~ r/ah. 



4.4 Thin disks with strongly magnetized coronae 

This is a variant of 



4.3 



in which most of the field lines are not open but re-attach to the disk at large 



distances Ar h (Galeev et al. 1978; Heyvaerts & Priest 1989). The vertical magnetic scale height is then 
H ~ Ar ^> h. The angular momentum flux carried through the corona, 



Jcoi-(r) 



dr 

h 



B r Bcf, 

An 



BrBtf, 
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can be larger than the flux within the gas layer by a factor ~ H/h, so that the effective value of a might be 
as large as r/h (for H ~ r) without fields exceeding equipartion. 

The evidence for magnetized coronae that dominate angular-momentum transport is suggestive but 
inconclusive. Local simulations of MR1 generally find that the scale height of the field exceeds that of the 
gas, but only by factors of order unity; they also find a ~ 10 -2 — 10 _1 rather than ~ r/h (Brandenburg et al 



1995 



Stone et al. 1996 



Miller & Stone 2000). Possibly, H/h is limited by the fact that the smallest dimension 



of the computational domain is < h. Global simulations of MRI have been performed for relatively thick 
disks only, so that it is difficult to distinguish scalings with h from scalings with r (Matsumoto & Shibata 



1997 



Hawley 2000 ) . Global simulations of thin disks in three dimensions may not be available for some time 



because of the very large numbers of grid cells needed to resolve both the disk and the corona. Merloni & 



Fabian (2001) argue that X-ray observations of accreting black holes (AGN and X-ray binaries) demand a 



strongly magnetized corona, at least in the innermost part of the disk. On the other hand, observations of 
eclipsing cataclysmic variables indicate that X-rays are emitted from the disk-star boundary layer rather 
than an extended corona (Mukai et al. 1997; Ramsay et al. 2001). 



4.5 Global spiral waves 

As is well known, a trailing m-armed spiral density wave 
E(r,0)=E o (r) + E m (r) cos (m(9 + julnr) 

exerts an outward (positive) gravitational angular-momentum flux 

r * n*GYLr> jjf 



(37) 



(Lynden-Bell & Kalnajs 1972). The above approximation is good for tightly- wrapped waves, fi ^> m ^ 0, 
which carry relatively little flux for a given density contrast E m /Eo. An exact formula for logarithmic spirals 
with E m oc r _3,/2 is 

r=-7t 2 GE 2 Jl r 3 m-^-K(n,m), (38) 



where K(p,m) is the Kalnajs function (Kalnajs 1971): 



r[i(m + i M +§)] 



r[|(m + i/x + §)] 



(real jj,). 



(39) 



With E m /Eo = 1, the largest ratio for which the surface density is everywhere positive, one finds that the 
torque is maximized at m = 1 and a pitch angle tan _1 (m//i) ~ 48^2, so that 



0.96l7tGEor 3 . 



(40) 



Suppose that the gravitational torque is balanced by the advection of angular momentum with the 
accreting gas. In other words, T « MQr 2 , so that there is no secular change in the angular momentum 
within radius r. The gravitationally-driven accretion speed is then 

r 

L max 



0.15Q~ 



(41) 



27tr 3 fiE 

in which Q has been calculated from the azimuthal average of the surface density, Eo. In principle therefore, 
accretion may occur at a significant fraction of the sound speed. But the existence of a selfconsistent wave- 
driven flow has been assumed rather than proved. Nonlinear single-armed spirals in keplerian disks have 
been found by Lee & Goodman (1999), but only for weak self-gravity (Q 1) and without dissipation or 



accretion. 

In steady accretion onto a central mass that dominates the rotation curve, the advected angular- 



momentum flux r = Mflr 



power-law index does not change the results (^o|)-(^l]) much. 



so E oc r 5//4 rather than r 3//2 . Presumably the slight change in 
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4.6 Clumpy disks 



Can QSO disks persist at Q <C 1 without fragmenting entirely into stars? This question is all the more 
urgent because of rather direct evidence for parsec-scale accretion disks in nearby AGN, if not QSOs, from 
VLBI observations of maser emission. If the nuclear disk of NGC 1068 is in a steady state, then the nuclear 
luminosity implies Q ~ 10 -3 at r ~ 1 pc (Kumar 1999). NGC 4258 is much less luminous, and estimates of 
M range from 7 x 10 _5 ao.iM© based on modeling the maser emission itself (Neufeld & Maloney 1995), to 
10" " 



Mq (Gammie et al. 1999; Kumar 1999) for an assumed central ADAF; at the former rate, Q 
outer edge of the masering region, ~ 0.2 pc (Maoz 1995), while in the latter, Q ~ 10 -2 . 



1 at the 



Kumar (1999) has suggested a clumpy rather than smooth disk, in which accretion occurs by gravita- 
tional scattering and physical collisions among clumps rather than an a viscosity. These clumps are supposed 
to be gas clouds rather than fully formed stars, in order to provide appropriate conditions for maser amplifi- 
cation. Although the model deals with the stability and accretion rate of the clumpy disk as a whole, it does 
not ask what prevents the individual clumps from collapsing. In the application to NGC 1068, the masses, 
radii, and surface temperatures of the clumps are quoted as M c ~ 1O 3 M0, R c ~ 0.1 pc, and T e s tC ~ 500 K; 
the virial temperature implied by this mass and radius is ~ 2000 K. These clumps would be moderately op- 
tically thick at the wavelength corresponding to T G ft, c . Their characteristic thermal time — the time required 
to radiate their binding energy — is 



tth,. 



GM?/R C 
4nRZaT* 



hence many orders of magnitude less than the orbital time (f2 _1 ~ 10 11 s at r ~ 1 pc) or interclump collision 
time. Of course the surface temperatures of these objects are assumed to be maintained by irradiation from 
the central source, but without a fast-responding feedback mechanism, the thermal equilibrium is unstable: 
if a clump starts to contract, its surface temperature will rise, and collapse will proceed on the clump's 
internal free-fall timescale (since this is larger than tth.c)- 

To answer the question raised above, we therefore believe that no QSO accretion disk, whether smooth 
or clumpy, can persist at Q <C 1. 



5 DISCUSSION 

Although the accretion velocity of a thin viscous disk is very strongly subsonic, with a Mach number 
M = v r /c s « ah/r, the discussion of §^ suggests that M might approach ~ 0.1 if accretion is driven by 
large-scale magnetic or gravitational fields rather than a local effective viscosity. Self-consistent solutions 
that achieve this bound would be very interesting to pursue. 

Even at near-sonic accretion speeds, QSO disks become selfgravitating at radii less than a parsec if 
there is no important source of heating other than dissipation of orbital energy. Combining eqs. (Q) & (^) 
with M = 2nrY,Mc B instead of the viscous relation (Q), and assuming that gas pressure dominates (/3 ~ 1), 
one finds that 

/ 7 \ 8/9 / \ -25/18 

-*B/l* n -l M 7/9-2/9 K 7/9 G -l ( _^B_) R -ll/9 f _T_} 

0.8 



1/4 mJ s \R S J 



-25/18 / M 7 5 ~2 X 1/9 



W*Rs J V llMl 1 

-25/18 f . .7 5 -2\ 1/9 



On the other hand, we have seen (§3.2) that when angular- momentum transport is viscous, then it is 
unlikely that stars can supply enough additional heat to stabilize the disk beyond one parsec, especially if 
viscosity is proportional to gas pressure rather than total pressure as viscous stability probably demands. In 
the latter case, even low-mass black holes embedded in the disk are probably inadequate. These statements 
assume that enough free gas remains in the disk to supply the central black hole at its Eddington rate. If all 
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the gas converts to stars, stability may result, but the quasar is quenched. Perhaps a combination of stellar 
(or embedded-black-hole) heating and a super-viscous accretion speed may allow an extended gravitationally 
stable disk; we hope to explore this possibility in a future paper. 

Given the serious theoretical difficulties of all proposed mechanisms for keeping Q > 1 at large radii, 
we are forced to take seriously the only remaining possibility: that QSO disks do not exist much beyond 
r sg ~ 10~ 2 pc [eq. (0)] — at least not in a state of centrifugal support, vertical hydrostatic equilibrium, and 
steady accretion. Yet the mass within this radius is smaller than that of the central black hole by a factor 
~ (h/r) ~ 10 -2 (since the midplane density « M/2nr 3 at Q — 1). Hence in order to grow the black 
hole, the disk must be replenished, either steadily or intermittently, by infall of low-angular-momentum 
material. In order that the gas not circularize outside 10 _2 r pc , its specific angular momentum must be 
< 70 Mg /2 k m s' 1 pc, or some three orders of magnitude smaller than that of most stars in ellipticals 
and bulges (Binney & Merrifield 1998). Such a small ratio may seem unlikely, but on the other hand, 
A^bh/A^buigc ~ 10 -3 (McLure & Dunlop 2001). So perhaps the QSO is fueled by the low-angular-momentum 
tail of gas that forms the bulge. This gas would arrive at the outer edge of the disk in a vertically broad 
infall, perhaps already carrying dust formed from metals injected by outflows from the bulge or QSO disk 
itself, and hence taking the place of a warped outer disk as the source of reprocessed infrared light. The 
picture this calls to mind is similar, except in scale, to the standard scenario for the formation of a protostar 
( |5hu et al. 1987[ ). 

I thank Charles Gammie, Pawan Kumar, Kristen Menou, Ramesh Narayan, and Bohdan Paczyhski for 
helpful discussions. This work was supported by the NASA Origins Program under grant NAG 5-8385. 
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APPENDIX A: VISCOUSLY HEATED DISKS 

For completeness, this appendix gives formulae for the midplane temperature (T), the surface density (S), 
and the gravitational stability parameter (Q) in a steady disk heated by viscous dissipation only, and cooled 
by radiative diffusion. 

Combining eqs. (jjj), and (^), and writing /3Cg = feT/m, where m ~ uih is the mean mass per gas 
particle, we have the radial dependence of £ & T: 

T = ( V /5 M 2 / 5 n 3 / s (Al) 
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g cm 2 . 



If viscosity scales with gas pressure (6 = 1) then eqs. ( |A2[ )-( |aT| ) do not depend on f3, which in any case is 
not an independent parameter: 

P _ Pgas _ 3efc B p_ _ 3c / fc B \ 1/2 ^i/ 2 S» 



1 - /3 Prad 4(Tm T 3 8(7 V ?Tt 7 ^ T 7 / 2 

this leads to 

fl (l/2) + (6-l)/10 / _ \ 21/20 

P ~ n aa _-l/10 4/5,-4/5 / r \ 



i-, - oa ^^'{wr s ) • ^ 

So the importance of gas pressure increases monotonically with radius in a steady, viscously heated disk. 
Using eq. (Al) to eliminate c a from eq. (^), 



i ^7/10^(76-125/10 ( _fcB \ G-l^-2/5^9/10 (A4) 



(47t) 3 / 



-2 7/10 fl (76-12)/10 / Sa±\ -3/10 ,,-13/10 / lO^Rs x ~ 
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